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Abstract 



CN ! The di-electron widths of ■0(4040), ■(/'(4160), and 0(4415), and their 

ratios are shown to be in good agreement with experiment, if in all 
' cases the S — D mixing with a large mixing angle 9 ~ 34° is taken. 

Arguments are presented why continuum states give small contribu- 
Qh| tions to the wave functions at the origin. We find that the y(4360) 

(— I ■ resonance, considered as a pure 3 -Di state, would have very small 

di-electron width, ree(i^(4360)) = 0.060 keV. On the contrary, for 
\ large mixing between the 4'^S'i and 2>^Di states with the mixing angle 

' 9 = 34.8°, ree('i/'(4415)) = 0.57 keV coincides with the experimental 

number, while a second physical resonance, probably y(4360), has also 
\ a rather large TeeiXi"^ 4400)) = 0.61 keV. For the higher resonance 

CN \ y(4660), considered as a pure 5^5i state, we predict the di-electron 

lO ' width ree(5^(4660)) = 0.70 keV, but it becomes significantly smaller, 

, namely 0.31 keV, if the mixing angle between the 5^51 and 4^Z)i 

(<~^ \ states 9 = 34°. The mass and di-electron width of the 6^5"! charmo- 

nium state are calculated. 



X. 

c5 ■ 1 Introduction 



Knowledge of the di-electron widths of higher charmonium states is important 
for many reasons. First of all, it can help to identify the nature of the newly 
discovered resonances with J^*-^ = 1 and distinguish between conventional 
cc mesons and, for example, tetraquarks which have much smaller di-electron 



^ Institute of Theoretical and Experimental Physics, Moscow, Russia. 
^Department of Physics and Astronomy, Vrije Universiteit, Amsterdam, The Nether- 
lands. 

■^Moscow Engineering Physics Institute, Moscow, Russia. 
* E-mail: badalian@itep.ru 



1 



widths pp. As shown in pp, di-electron widths of compact four-quark systems, 
hke ccgg, with J^^ = 1 appear to be two orders smaller than those of 
conventional cc mesons. 

For higher charmonium states, which lie above the DD* (or D*D*) thresh- 
old, their wave functions (w.f.) can be strongly affected by the coupled- 
channel and threshold effects, being in general very complicate functions. 
Via open channels the w.f. of a resonance contains admixtures of other 
states with the same quantum numbers and also a contribution from contin- 
uum state(s). To define such a w.f. one needs to formulate the relativistic 
many-channel Hamiltonian in QCD, even if in some approximation [2J. In the 
simplest case, the continuum part consists of two open-charm mesons and to 
some extent this continuum part can be considered as a particular case of a 
four-quark system, ccqq. Since the contribution of any four-quark state to the 
w.f. at the origin is much smaller than that of a meson, it can be neglected. 
This effect occurs because the probability to collect four (and even three) 
particles at the origin is much smaller than for two particles. This fact does 
not exclude that at larger distances the continuum part can give an essential 
contribution to the w.f. and even dominates asymptotically. Just such a 
coupling to the continuum provides a shift down of the mass of the P-wave 
heavy-light mesons [2]. Therefore the w.f. at the origin of a higher vector 
resonance, which we are interested in here, can be calculated taking into ac- 
count only the mixing between those vector states which masses, defined in 
single-channel approximation, have close values. Study of the charmonium 
spectrum shows that the (n + 1) ^5*1 and n ^Di states (n > 3) have small mass 
differences, < 60 MeV, which decrease for higher radial excitations. On the 
other hand the mass differences between neighbouring ^Si states is of the 
order of several hundreds of MeV, so in first approximation mixing between 
these states can be neglected. Such a representation of the w.f. at the origin 
can be tested via concrete predictions for the di-electron widths of different 
vector states in heavy quarkonia. 

In this picture the S — D mixing between higher resonances can be con- 
sidered in the same way as it has been done for iIj'{3686) and ip"{3770) [5], 
1^, where the mixing angle 6 = (12 ± 2)° is extracted from the ratio of their 
di-electron widths. Here we show that for higher vector states the S-D mixing 
is significantly larger and the mixing angle 6 ~ 34°. Just for such an angle 
the di-electron widths of '?/'(4040) and ijj{4:lQ0) turn out to be almost equal, 
as in experiment [5] -[9]. 

There are also other arguments in favor of a large S — D mixing. It is 
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known that the experimental di-electron widths of '?/'(4040), '?/'(4160), and 
ip{4415) are significantly smaller than the potential model predictions in 
single-channel approximation [10]- [13]. Moreover, even with the use of a 
many-channel w.f. at the origin (calculated in [12] for the Cornell coupled- 
channel model [H]) the di-electron widths of higher charmonium states ap- 
peared to be considerably larger than in experiment. It is well-known that the 
QCD-motivated gluon-exchange (GE) potentials (with the strong coupling, 
possessing the asymptotic freedom property) have much smaller w.f. at the 
origin than those with ag = const [13]; nevertheless, even such potentials give 
di-electron widths of ip{404:0) and ip{4415) which are still 50 — 70% larger 
than the experimental numbers [in]-[I2]- In the recent paper [TJ] the y(4660) 
resonance [IB], considered as a 5^ Si state, has also large ree(5^S'i) = 1.34 
keV, which is even significantly larger than Ted'ip (4415)) = 0.58(7) keV. 

We assume here that the relatively small values of the di-electron widths 
of the n ^Si states and the rather large widths of the n ^Di states, initially 
considered as pure states, occur mostly due to S — D mixing. Here we do not 
use the same assumption as made in [TT], where to reach agreement with ex- 
perimental widths, "total" screening of the GE interaction at large distances 
has been supposed: such an assumption cannot explain why ree(V'(4160)) is 
large, and also has no deep theoretical grounds. 

Three experimental facts point to a possibly large mixing between the 
{n + 1) ^Si and n ^Di vector charmonium states: 

1. The measured di-electron width of ip{4160), which is usually considered 
as the 2^Di state, is large [8], [9]: 



Namely, it is only 5 — 10% smaller than the di-electron width of iIj{4040) 
and ~14 times larger than the width calculated here for a pure 2^Di 
state: Teei^^Di) = 0.061 keV (see Section 3). It is also about three 
times larger than T ^ei^p" (3770)) = 0.248(6) keV. 

2. On the contrary, the experimental width of ?/'(4040), considered usually 
as the 3 ^Si state. 



appears to be almost two times smaller than for a pure 3 ^Si state [12] . 



ree(V'(4160)) = 0.83 ± 0.06 keV. 



(1) 



Fee (^A (4040)) = 0.86 ± 0.07 keV, 



(2) 
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This situation can be resolved if the S-D mixing between these two 
states is taken into account. For levels above the D*D* threshold such 
a mixing can occur owing to short-range tensor forces (it gives a rather 
small effect) and the influence of open channel(s). Since at present 
there is no dynamical calculation oi S — D mixing, the influence of open 
channels can be taken into account in a phenomenological way through 
the introduction of a mixing angle, as for '?/''(3686) and '?/'"(3770) [3], 

3. The third fact refers to the di-electron width of '?/'(4415). If this res- 
onance is considered as the 4^S'i state, then potential models give di- 
electron widths in the range 1.1 — 1.5 keV [10], [12], [13], which are 
almost two times larger than in experiment [8], [9]: 

ree(^(4415)) = 0.58 ± 0.07 keV. (3) 

Such a decrease of the di-electron width could occur via mixing with 
a still unidentified 3 ^Di state, which in single-channel calculations has 
mass M{?,D) = 4.470(10) MeV, while M{4S) = 4420(10) MeV, i.e., 
these two masses are rather close to the masses of the physical reso- 
nances '?/'(4415) and y(4360). One may expect that these 4^5*1 and 
3 ^Di states could be strongly coupled to the 5- wave decay channels, 
like Di(2420)D*(2010), L'*(2400)L'*(2010), and L'*o(2317)L':(2112), 
and due to this coupling the 4 ^5*1 and 3 ^Di levels are mixed and acquire 
hadronic downward mass shifts, which are typically ~ 40 — 60 MeV [H]. 
As a result, one of the shifted physical states goes over into the conven- 
tional '?/'(4415) charmonium, while the other one can possibly be identi- 
fied with the y(4360) resonance, recently discovered by the Belle Col- 
laboration [TB]. (In our analysis here, the Belle resonance y(4360) with 
r = 48(15) MeV [T6] and the wide resonance F(4324) with T = 172 
MeV, observed by the BaBar collaboration [T7j, are considered to be the 
same). Then the di-electron width of '?/'(4415) is calculated here, taking 
into account large S — D mixing, while for the analysis it is inessential 
from which state, 4^5*1 or 3^Di, the resonance '?/;(4415) originates. We 
show that for 6 = 34° the di-electron widths of both physical resonances 
have close values: ree(l"(4360)) ~ ree(V^(4415)) = 0.58 keV. 

In our picture it is convenient to define the mixing angle between higher 
vector states from the ratio of the di-electron widths, as in [5], [1]: in this 
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case the QCD factor occuring due to radiative corrections (see Sect.H]), is 
cancelled in the ratio. From such an analysis a large mixing angle is extracted 
and the absolute values of ree(V'(4040)) and Fee (^/^ (4160)) are obtained in 
good agreement with experiment if the same QCD factor jSy = 0.63 is taken 
for all higher states. 

Notice that the mass of the 5=^^i state, M{5^Si) = 4640(10) MeV, has 
been predicted in [12], before the Belle resonance F(4660) was discovered 
[T6]. This resonance and its radiative transitions were studied in detail in [T5] 
giving Fee (y (4660)) = 1.34 keV. In our calculations Fee (>" (4660)) strongly 
depends on a possible admixture of the 4 ^Di state and is considerably smaller 
than in [15]: F(y(4660)) = 0.70 keV, if this resonance is a pure 5^5*1 state, 
and about two times smaller, Fee(y(4660)) = 0.31 keV, if the mixing angle 
between 5 ^5*1 and the unobserved 4 "^Di state (with mass ~ 4700 MeV) is 
34°, the same as for '?/'(4415). 

2 The masses of the J^*" = 1 charmonium 
states 

The hyperfine (HF) and fine-structure splittings of higher radial excitations 
are small (< 20 MeV) [TT], [TB], therefore their masses practically coincide 
with the centroid masses, McogijiL), which we need to determine with good 
accuracy. To calculate them we use here the relativistic string (RS) Hamil- 
tonian with universal (for all mesons) interaction [19] , [20] . For charmonium 
one contribution to the mass formula, namely the small string correction 
(< 5 MeV ) for the states with L 7^ can be neglected, while the self-energy 
correction, ~ —20 MeV, is taken into account here. 

In heavy quarkonia the mass Mcog{nL) is just given by the eigenvalue 
(e.v.) of the spinless Salpeter equation (SSE) [20] : 



Here we also use the RS Hamiltonian written in the Einbein approximation 
(EA) [21] . In this case the spin-averaged mass can be presented as: 




(4) 



Mcog(nL) = UnL 



+ 




(5) 



5 



where the e.v. Eni are the solutions of the so-called einbein equation 

" p2 1 

\- Vo{r) iPnL{r) = EnL'f nL, (6) 

which together with the mass uJnL should be defined in a selfconsistent way: 

For the n^Di state the small HF contribution to the mass will be ne- 
glected here, i.e., its mass M{n^Di) = Mcog(nD), while for the n^Si states 
we still keep the small HF correction: M{n^Si) = Mcog{nS) + j6uF{nS) with 
5uF{nS) = M{n^Si) - M{{n^So). The values of dnF{nS) = 48(48), 16(20), 
12(16), 6(10) MeV (n = 2, 3, 4, 5), calculated in [H] and p] (in parentheses), 
are used here. 

Our calculations are performed with the universal potential VB{r) from 
EOl, [22]: 



VB(r) = a(r)-r-^^, (8) 

where the vector coupling asir) is taken in two- loop approximation: it has 
the asymptotic freedom behavior at small r, freezes (saturates) at large r, 
and depends on the number of fiavors Uf. For charmonium we use nj = 4 and 
the QCD constant = 254 MeV, which gives the vector QCD constant 
Ay(nf = 4:) = 1.4238 ■ A^sinf = 4) = 360 MeV [22]. The freezing (critical) 
value of aB{r) is expressed through Ay and the so-called background mass 
Mb = 1.0 GeV: 

aB{r ~^oo)= aB{q = 0) = -— 1 - ^— ^ = 0.546, (9) 







where to = In (^) 2, (3o = 11 - |n/, and (3i = 102 - fuf. 

For low-lying states which have relatively small sizes (with r.m.s. radius 
R{nL) < 0.8 fm), a linear confining potential with constant string tension, 
cr = (To = 0.18 GeV^, can be used [23]. However, for higher states, which lie 
above open thresholds and have large radii R{nL) > 1.0 fm, it is important to 
take into account the creation of virtual light-quark pairs (gg), even in single- 
channel approximation. Due to virtual loops the surface inside the Wilson 
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Table 1: The charmonium masses M{n^Si) (in MeV) for the potential (JH]). 



state 






BGS[18] 


exp.[9j| 




rric = 1.425 GeV 


me = 1.410 GeV 


TTic = 1.4794 GeV 




IS 


3105 


3095 


3090 


3097 


2S 


3678 


3682 


3672 


3686 


3S 


4078 


4096 


4072 


4039(1) 


AS 


4398 


4426 


4406 


4421(4) 
4361(18)^) 


5S 


4642 


4672 




4664(16)^^ 


6S 


4804 


4828 







" The self-energy corrections to the masses, AsE{nS) ~ —20 MeV, are taken 
into account both in the relativistic case (SSE) and in the einbein approxi- 
mation. 
Belle data [16] 



loop decreases, making the string tension dependent on the QQ separation r 

a(r) = ao(l-7/(r)). (10) 

Such a flattening of the confining potential is common to all mesons of large 
sizes, and therefore for charmonium the form and parameters of such amod- 
ified string tension can be taken from the analysis of the radial Regge tra- 
jectories for light mesons |24j : 

7 = 0.40; /(r^0)=0, /(r ^ oo) = 1.0. (11) 

As shown in |23], due to flattening the masses of all higher levels are shifted 
down and these mass shifts increase with n. For example, the shift of the 
5^5*1 state reaches ~ 100 MeV [23j. In Tables [T] and [2] we give the masses of 
pure n^Si and n^Di states for the potential ([HD, which are calculated using 
the SSE and the EA ([H]), and also the masses calculated in [T8], where in 
the Cornell potential a constant coupling is used, equal to our freezing value 
(Q. From Table [1] one can see that in our calculations the mass of the 5^5*1 
state agrees with that of the y(4660) resonance [16]. For the 6^5*1 level the 
predicted mass is M(6^5i) = 4815(15) MeV. 
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Table 2: The charmonium masses M{n^Di) (in MeV) for the potential 



state 






BGS[18] 


Exp. [S] 




rric = 1.425 GeV 


me = 1.410 GeV 


TTic = 1.4794 GeV 




ID 


3800 


3779 


3806 


3770(3) 


2D 


4156 


4165 


4167 


4159(3) 


3D 


4464 


4477 




4421 

4361''), 4324^) 


AD 


4690 


4707 






5D 


4840 


4855 







" See the footnote " to Table □ 
^ See the footnote to Table [T] 
" BaBar data [HI 



As seen in Tables [T] and [21 the masses M((n + 1)5*) and M{nD) are close 
to each other, even in single-channel approximation. The difference between 
them decreases for larger radial excitations, so that for = 5 it is only ~ 30 
MeV. It is worthwhile to notice that besides the "correlated" mass shifts 
of higher levels — due to virtual pair creation — some levels, which lie near 
thresholds, can have additional downward shifts due to strong coupling to 
a continuum channel. We denote these mass shifts as decay-channel (DC) 
shifts; they can be calculated only within a multi-channel approach. The 
masses of the 25", 45", and 2D levels calculated here (see Tables [T] and [2]) 
agree with experiment within 20 MeV, i.e., they have essentially no DC 
shifts. However, the mass of the 3 ^5*1 level is ~ 40 MeV larger than the 
experimental one, because this level can be affected by the D*D* channel 
and we estimate its DC shift as ~ 40 MeV. Therefore, the masses predicted 
in our paper, have different accuracies, which is better than 20 MeV for the 
levels without DC shifts, and than 40 MeV for the levels strongly coupled to 
nearby continuum decay channels. 

Thus from our analysis we conclude that the mass difference, 

Ar,M = M{n^Di) - M{{n+l)^Si), (12) 

decreases from the value A2M(exp) = 120 MeV for = 2 to ~ 30 MeV for 
n = 4. Therefore higher levels are almost degenerate and the S — D mixing 
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for them, as well as the mixing angle, become larger than for ip' and ip" . 
Also in single-channel approximation Mcog{{n +1)5') is always smaller than 
Mcog(nD). 



3 Vector Decay Constants 

The decay constants of vector (V) mesons are calculated here using the an- 
aljd;ical expressions derived in [25] in the framework of the Field Correlator 
Method, where fv{nS) and fv{nD) are given by 

,J.(.„S).12^ft.^i§i(f (13) 

My{nS) 71 Mv{nS) 

where the relativistic factor is defined as 

iv[nL) = ^-2-^ . (14) 

The same expression ( IT3il can be applied to the decay constants fv{n^Di), 
if in (fT3|) the w.f. Rnoi^) is defined as in (fT5l) below. 

All matrix elements (m.e.) which are needed to calculate the decay con- 
stants for the nS and nD states, are given in Tables [3] and HI An interesting 
fact is that for the (n + 1)5 and nD states the m.e. like UnL, (p^)nL, and 
(^y(?7.L) coincide with an accuracy better than 1% and therefore the difference 
between fv{{n + 1)5) and fv{nD) comes only from their w.f. at the origin 
and the small differences in their masses My{nL). 



Table 3: The matrix elements Uns (GeV), (p^) ((GeV/c)^), and the w.f. at 
the origin Rns{0) (GeV^^^) (no mixing) for the potential ([8]) 



State 


cOnS (GeV) 


Rns{0) 






1^ 


1.59 


0.905 


0.541 


0.929 


2S 


1.65 


0.767 


0.722 


0.910 


3S 


1.69 


0.714 


0.882 


0.899 


AS 


1.71 


0.655 


0.947 


0.894 


55 


1.66 


0.531 


0.775 


0.908 


65 


1.63 


0.445 


0.665 


0.916 
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Table 4: The matrix elements UnD (GeV), (p^) ((GeV/c)^), the w.f. at the 
origin RnD{0) (GeV^/^), and the second derivative KoiO) (GeV^/^) for the 
potential ([8]). 



State 


LUnD (GeV) 


<d(0) 


RnD{Q) 






ID 


1.65 


0.145 


0.095 


0.721 


0.909 


2D 


1.69 


0.213 


0.132 


0.881 


0.899 


3D 


1.71 


0.248 


0.150 


0.939 


0.893 


AD 


1.65 


0.221 


0.144 


0.745 


0.911 


5D 


1.64 


0.206 


0.135 


0.682 


0.912 



The w.f. at the origin i?„£)(0) is defined here as in [26], being expressed 
via the second derivative i?„^(0): 

Rnom = J R'UO). (15) 

It is interesting that Rns{0) and i?nD(0) have different behavior for growing 
n: while Rns{0) decreases for higher radial excitations, the second derivative 
R'^^{0) and RnoiO) grow with increasing n, if a linear confining potential is 
used. For the flattening potential used here, Rns{0) decreases even faster, 
while Rnoi^) increases for the 25" and 3S states and then practically saturates 
for higher levels. This growth of i?„£)(0) is possibly one of the reasons why 
higher radial excitations have large S — D mixing. 

In Table [5]the decay constants of the charmonium states with J^^ = 1 
are given in two cases: without and with S — D mixing. The mixing angle 
between the 2 ^5*1 and 1 ^Di levels has already been calculated in [3] and [1] , 
and also in [23], where the mixing angle ^ = 11° has been extracted from the 
ratio of experimental di-electron widths, ree(V''(3686)) and Tee{ip (3770)). 
Other mixing angles are calculated below. The matrix elements and other 
numbers in Tables |3] and H] are calculated here with the use of the EA equation 
(ED, which provides regular behavior of the w.f. at the origin, in contrast to 
the SSE for which the S"— wave w.f. diverges at the origin and needs to be 
regularized. 

As seen in Table |5] the decay constants fv{nD) are very small, ~ 50 MeV 
{6 = 0) for pure D-wave states. However, if S ~ D mixing is large {6 ~ 34°) 
the decay constants fv{(^) (^ > 2) of physical, "mixed" states appear to be 
practically equal for '?/'(4040) and iIj{41Q0), and also for ?/^(4415) and the still 
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Table 5: The decay constants fv{n^Si) (in MeV) without and with S — D 
mixing. 



9 = 




n = 1 


n = 2 


n = 3 


/v((n+l)35i) 


373 


329 


288 




45 


60 


66 




6 


11° 


34.8° 


34° 


fvi^s) 


357 


236 


202 




115 


234 


217 



unidentified second charmonium state, which originates from the 2D level 
and it is denoted below as ip^Mld) (although its mass is ~ 4470(10) MeV in 
single-channel approximation). For large S — D mixing all decay constants 
lie in the range 220 ± 20 MeV {n = 2,3). Precisely this fact provides close 
values of the di-electron widths of the (n + 1)S and nD states. 



4 Di-electron widths 

In [23] the di-electron widths: Tee{J/i'), Tedi'iSGSG)), and Tef>{ip{3770)) have 
been calculated with high precision, < 5%, using the theoretical formula 
where the di-electron width is expressed via the decay constant (IT^ . con- 
taining the relativistic correction C,v, and includes QCD radiative corrections 
(this expression is the relativistic generalization of the van Royen-Weisskopf 
formula [27] in the framework of Field Correlator Method). The QCD cor- 
rection, known in one-loop approximation, enters as the multiplicative factor 
denoted here as /3y = 1 — |^as(My). Then 



reein'S,) = ^^flsPv = -^\RnSm%sPv, (16) 



4:7ieW „o . 4e^a^ 



reein'D,) = ^--f^f^^Py = -^\R^^(0)\%j,f3y. (17) 



The w.f. at the origin Rnoi'^) in (jlTl) has been defined in (fTSl) and the average 



kinetic energy uOni = (yP^ + "mDnL, calculated from equation ([7j), plays the 
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role of a constituent quark mass being different for different nL states. The 
uJnL have the following characteristic feature: For a linear confining potential 
with a = (Jq = const, it grows for higher radial excitations, while for the 
flattening potential first it grows for n = 2, 3, and 4 and then saturates 
around the value Uni ~ 1-65 GeV for n > 5 (the values Uns, ^no are given 
in Tables E] and SD. 

The expressions ffT^ and flTTl) contain the c-quark charge Cc = 2/3, a = 
1/137, and the mass MnsiMno) of the n^Si{n^Di) vector mesons. The w.f. 
at the origin Rns{0), RnviO), and i?^^(0) are given in Tables [3] and HI 

As we discussed in the Introduction and Section 2, one may expect that 
the physical ■i/'-mesons represent a mixing of the {n + 1)S and nD states with 
close mass values, and our goal here is to determine the mixing angle between 
higher radial excitations. To this end we introduce the w.f. of the physical 
■j/'-mesons, denoted here by v'nS'(O) and (pnoi^), where the symbols nS and 
nD simply remind about the origin of those states: 

<PnsiO) = COS 9nRns{0) - sin 9nR(n-l)D{0), 

VnoiO) = COS 9nR(n+l)S (0) - sin 9nRnD{0). (18) 

The di-electron widths of the ■j/'-mesons are expressed via the physical w.f. at 
the origin (fT8|) in the same way as in (fT6|) and (ITTl) . For a given mixing angle 
9 the w.f. at the origin (fTS!) are easily calculated through the w.f. Rns{0) 
and i?„£)(0) for pure S- and D-wave states (they are given in Tables [3] and 
S}. 

Table 6: The wave functions at the origin Lp(^n+i)s{0) and ifnDiO) in GeV'^/^ 
of the physical states for ri = 1, 2, 3, 4". 



n 


1 


2 


3 


4 


9 


11° 


34.8° 


34° 


34° 


(P{n+l)s{0) 


0.735 


0.511 


0.459 


0.360 




0.240 


0.516 


0.491 


0.416 



" The uncertainty in the mass value used gives rise to a theoretical error less 
than 1%. 
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5 S^^i - 2 3i:)i mixing 



To determine the mixing angle between the 3^5*1 and 2^Di states we use 
here the ratio of the di-electron widths of the tp{AOAO) and iplAlQO) mesons 
as in [4J and [22]. This ratio does not depend on the QCD factor jSy and 
the experimental di-electron widths are given in ([1]) and ([2]) with their ratio 
close to unity. Such a large ratio turns out to be possible only if the mixing 
angle between 3 ^Si and 2 ^Di states is large. 

Taking in the ratio the w.f. at the origin f[T5]) . which are expressed via the 
numbers Rns{0) and i?„D(0) from Tables [3] and IH one can extract the mixing 
angle 62 = 34.8° and determine the physical w.f. at the origin, as well as the 
decay constants of ip{4040) and iIj{4160). Then for both charmonium states 
the di-electron widths appear to be in precise agreement with experiment 
(see (l2Til ). if the QCD factor Pv = 0.63 is taken. This value of Pv is smaller 
(i.e. the radiative corrections are larger) than for J/ip, '?/'(3686), and ■ip{3770), 
where in all cases the larger value /3y = 0.72 gives precise agreement with 
experiment [23]. 

The mixing angle between the {n + 1) ^5*1 and n^Di states, denoted here 
as 6n, can be calculated if at least one of the di-electron widths is known from 
experiment. For the 35* and 2D states both di-electron widths are known 
and 62 is easily determined. It is important to notice that for a pure 2 ^Di 
state the di-electron width is very small: ree(2 ^Z}i)=0.059 keV, i.e., ~ 14 
times smaller than the experimental number ([1]), and one can expect large 
mixing between the 3 ^Si and 2 ^Di states. Such a large mixing can occur 
via the nearby open D*D* channel and partly through short-ranged tensor 
forces which, however, do not provide a large mixing angle, ^(tensor) < 7°. 
On the contrary, from the ratio: 

r„W.(4040))^ (19) 
' r„(i/.(4160)) ^ ' 

one obtains two solutions with a large magnitude of 6'2: a positive and a 
negative one: 

02 = 34.8° or 02 = -55.7°. (20) 

For these angles and using (fT8l) the physical w.f. (f{ip{4:04:0),r = 0) = 0.511 
GeV3/2 and ¥^(^^(4160), r = 0) = 0.516 GeV^/^, appear to be almost equal. 
Then from (ITBl) and (ITTI) with (3v = 0.63 one calculates the following di- 
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electron widths: 



ree(^(4040)) = 0.87 keV, ree(V'(4160)) = 0.83 keV, (21) 

which just coincide with the central values of the experimental values, ([1]) 
and ([2]). The QCD factor P = 0.63 extracted simultaneously, corresponds to 
the strong coupling ^^(Mv) = 0.217. Later this value of /3y = 0.63 is used to 
determine 6n for higher excitations {n = 3,4). Notice that the same mixing 
angle 62 = 35° (or 62 = —55°) has been obtained in the analysis of '?/'(4040) 
and V^(4160) [28]. 



6 Large mixing between and 3^Di states 

In constituent quark models (in single-channel approximation) two vector 
states, 4^5*1 and 3^Di, are expected in the mass region around 4.4 GeV 
(see Tables [1] and [2]). Our calculations give the masses M(4^S'i) ~ 4.42 
GeV and M{3^Di) ~ 4.47(1) GeV with their mass difference ~ 50 MeV. 
However, one cannot exclude that due to strong coupling to the D*Di{24:20) 
and D*D;{2AQ0) channels the Si and 3 3L'i states are mixed and the mass 
of one or probably both states is shifted down. Then one of these mixed 
(physical) states can be identified with ■?/'(4415) and the other one with the 
newly discovered resonance y(4360) [16j. (Note that in charmonium the 
values of the DC shifts are typically ~ 40 MeV [Ti]). 

From experiment only the di-electron width ree(^/'(4415)) is presently 
known, ree(V'(4415)) = 0.58 ±0.07 keV, while for the F(4360) resonance two 
possible numbers have been measured for the product [16], 

B(r(4360) - ^(2^)vr+vr-) x 1^^(^(4360)) = Jjj JJ;^ ^ (22) 

Still, this restricted information allows one to draw an important conclusion. 
First, for pure 45" and 3D states (6^3 = 0) with the w.f. at the origin R4,s{0) = 

0. 65 GeV^/^ and RsniO) = 0.150 GeV^/^ (from Tables [3] and H, their di- 
electron widths are the following: 

ree(4=^^i) = 1.19 keV, ree(3^Di) = 0.06 keV, (^ = 0). (23) 

1. e., ree(45') is two times larger than the experimental number ([3]) while 
ree{3D) is small. To reach agreement with experiment for ?/'(4415) we need 
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to take a large mixing angle, namely 63 = 34°, as for the 3S — 2D mixing, 
for which 

theory 

0.57 keV (24) 

is completely in agreement with the central experimental value ([3]). 

Then for the same angle the di-electron width of the second physical state, 
which can be denoted as ?/'(4470) (in many-channel approximation its mass 
may be smaller), appears to be ten times larger than for a pure 3^Di state: 

ree(V^(4470)) = 0.63 keV. (25) 

Moreover this width is even slightly larger than that of '?/'(4415) (here we 
take j3v = 0.63 as for the 35* and 2D states). Since in this case the di- 
electron widths coincide within 10% accuracy, in the framework of the single- 
channel approximation it is difficult to decide which of these states should 
be identified with ■?/;(4415) or with y(4360). From the experimental value 
f l22|) and the di-electron width fl28l) one obtains an estimate of the branching 
fi(r(4360) ^ ij{2S)n+n-), 

S(F(4360) ^ V(25)7r+7r-) ^ (1.6 ± 0.6)%, (26) 

which is rather large. Thus for large mixing angle (^3 = 34°) one cannot 
a priori conclude which resonance, 'i/'(4415) or y(4360), originates from the 
4^5"!. A decisive test for their identification could come from the study of 
their radiative transitions to the Xcj states, because radiative transitions 
are very sensitive to the mixing angle, as has been shown for iIj'{3686) and 
?/^"(3770) in m, 

7 y(4660),y(4815) 

The higher level y(4660) with M = 4660 ± 16 MeV, recently discovered in 
|16j . has a surprisingly small width, F = 48 ± 18 MeV. This state lies close 
to the 5- wave threshold D*D si{2535) (with the threshold mass Mth = 4647 
MeV) and to the P-wave threshold D{2^Si)D* with Mth = 4647 MeV (our 
calculations give the mass M{D{2^Si)) ~ 2640 MeV). Our predictions for 
the masses of the 5 ^5*1 and 4^Di states in single-channel approximation (see 
Tables [T] and [2]) are 

M{5^Si) = 4655(15) MeV, M{A^Di) = 4700(10) MeV, (27) 
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where the theoretical uncertainty is taken into account. These masses differ 
only ~ 50 MeV and large mixing between the two states can be expected. 
Unfortunately, at present their di-electron widths remain unknown, and here 
we calculate them, assuming that the QCD factor /5y = 0.63 and 6^ = 34° as 
it takes place for 45" and 3D states, and also first consider pure 55* and AD 
states {04 = 0), for which 

Teei^^Si) = 0.73 keV, Tee{A^Di) = 0.055 keV, {64 = 0) (28) 

i.e., ree(5S') is rather large, being even larger than ree(4S') ([3]), On the con- 
trary TeeiAD) is very small. 

For large S — D mixing with 64 = 34° (like 62 and ^3) we obtain 

ree(V^(4660)) = 0.32 keV, {64 = 34°) (29) 

which is two times smaller than fl28l) . For the second state, denoted as 
7/i(4690) we find. 

Fee (^/i (4690)) = 0.45 keV, (^4 = 34°) (30) 

the width appears to be eight times larger than for the pure AD state in (128!) 
and even larger than Fee(^/'(4660)). Notice that equal widths are obtained 
for a bit smaller angle, 64 = 30°: 

Fee(V'(4660)) = Fee(V'(4690)) = 0.39 keV. (31) 

We predict also the 6 ^Si state although this state has very large r.m.s. radius, 
R = 2.5 fm, even in closed-channel approximation. Its mass is M{6^Si) = 
4815 ± 15 MeV and T^^{Q^Si) = 0.20 keV for 65 = 34°. The existence of so 
high a resonance would be important for the theory. 

8 Conclusions 

We have studied the di-electron widths of higher n'^Si and n^Di radial ex- 
citations in charmonium and shown that 

1. The almost equal values of Fee(4040) and Fee(4160), as well as the 
small value of Fee (4415), can be explained, if large S-D mixing between 
{n + 1) ^Si and n ^Di states takes place. 
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2. For ■0(4040) and ■0(4160) precise agreement with experiment is obtained 
taking the mixing angle 82 — 34.8°. 

3. For ■0(4415) the calculated di-electron width coincides with the central 
value of the experimental width for a mixing angle 9^ — 34°. 

4. In all cases the QCD radiative corrections appear to be important and 
the same strong couphng as{Mv) — 0.217 is taken, giving ~ 30% effect. 

5. In the single-channel approximation used here DC mass shifts (due to 
strong coupling to a nearby threshold) cannot be calculated. Therefore 
it remains unclear which physical resonance, ■0(4415) or the Belle res- 
onance y(4360), corresponds to the 3^Di(4^S'i) state. For both states 
we predict close values of their di-clcctron widths: ree('0(4415)) = 0.57 
keV and ree(5^(4360)) = 0.63(7) keV (they coincide within the experi- 
mental error). 

6. Assuming that the 5^5*1 and A^Di states have also large S — D mix- 
ing, with = (32 ± 2)°), we obtain: ree(V^(4660)) = 0.35(4) keV, 
ree(V^(4690)) = 0.40(5) keV. 

7. One cannot exclude that a 6 ^Si state also exists, for which wc predict 
the mass M{6S) = 4815(20) MeV and di-electron width Fee = 0.20 
keV. 
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